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We discuss the importance of the real part of the Feynman- Vernon influence action for the analysis 
of dephasing and decay near the ground state of a system which is coupled to a bath. Using 
exactly solvable linear quantum dissipative systems, it is shown how the effects of the real and 
the imaginary part (describing friction and fluctuations, respectively) may cancel beyond lowest- 
order perturbation theory. The resulting picture is extended to a qualitative discussion of nonlinear 
systems and dephasing of degenerate fermions. We explain why dephasing rates will, in general, 
come out finite at zero temperature if they are deduced from the imaginary part of the action alone, 
a procedure which is reliable only for highly excited states. 



I. INTRODUCTION 

The Feynman- Vernon influence functionalQi plays a 
prominent role in discussions of dephasing that aim to go 
beyond a simple master equation approach. In principle, 
the influence functional takes into account every effect of 
the environment (the "bath") on the system under consider- 
ation. This includes heating, dephasing, friction and renor- 
malization effects (changing the external potential or the 
effective mass of a particle) . Its popularity arises not only 
from the fact that it constitutes an exact approach, but also 
from the direct physical meaning which it acquires in some 
situations. In particular, this holds for typical interference 
experiments, where two wave packets describing a single 
particle follow two different trajectories in order to be re- 
combined later on. For such a case, the influence functional 
is equal to the overlap between the two different bath states 
which result due to the particle moving along either one of 
the semiclassical trajectories. If the bath can distinguish 
between the two paths, it acts as a which-way detector, and 
the diminished magnitude of the overlap directly gives the 
ensuing decrease of the visibility of the interference pat- 
tern. Writing the overlap in the form of an exponential, 
exp(iSi?, — Si), one can obviously conclude that only the 
imaginary part Si of the "influence action" is responsible for 
this decrease and, therefore, Si alone describes dephasing 
in such a situation. Furthermore, if the system is subject 
to a fluctuating external classical field, only Si remains, 
while Sr vanishes identically. This is because Si is due 
to the fluctuations of the bath, while Sr stems from the 
back-action of the bath onto the system (including friction 
effects), which is absent for an external noise field. The ap- 
proach of evaluating Si along semiclassical trajectories has 
been successfully applied .-to the calculation of dephasing 
rates in many situationsBoB. 

However, near the ground state of a system, an analysis 
along these lines is likely to fail. Qualitatively, this may 
already be deduced from the fact that dropping Sr means 
replacing the environment by an artificial fluctuating classi- 
cal field whose correlation function includes the zero-point 
fluctuations of the original quantum bath. Therefore, even 
at zero temperature, this fluctuating field will, in general, 
heat the system, which directly leads to dephasing. The 
role of Sr is to counteract this effect. In this article, it is 
our aim to display explicitly, in a detailed manner, the ne- 



cessity of keeping the real part Sr of the influence action in 
discussions of dephasing and decay near the ground state 
of a system. 

This issue derives its importance partly from the fact 
that the evaluation of Si along semiclassical trajectories 
has proven to be an efficient way of extragtipg dephasing 
rates in the problem of weak localization&a in the limit 
of high temperatures, when zero-point fluctuations of the 
bath may be omitted. In contrast, at low temperatures, 
a single-particle semiclassical calculation may become in- 
valid, since it neglects the Pauli principle which is known 
to play an important role for the inelastic scattering of elec- 
trons and is not included in the Feynman- Vernon influence 
functional. However, recently an extension of the influence 
functional to the case of a many-fexmipn system has been 
derived, using an exact procedureSCIS and including the 
Pauli principle. This permitted a discussion of dephasing 
in a disordered metal even for the case of low tempera- 
tures. Following the general strategy of earlier workso deal- 
ing with the high-temperature case, the dephasing rate was 
deduced from Si in a semiclassical calculation and found 
to be finite even at zero temperature. Since the "orthodox" 
theoryQ had predicted a vanishing rate in the limit|X|=j|-Q, 
the new results prompted considerable criticismE3liM3llj, 
which mostly emphasized technical aspects of impurity av- 
eraging or used perturbation theory to arrive at different 
conclusions. In the present work, we want to clarify some 
essential aspects of the roles of Sr and Si, using physically 
much more transparent exactly solvable models. 

We want to emphasize that "zero-temperature dephas- 
ing" as such is perfectly possible: If one prepares a system 
in any superposition of excited states and couples it to a 
bath, it will, in general, decay towards its ground state by 
spontaneous emission of energy into the bath (at T = 0). 
This will destroy any coherent superposition, thus leading 
to dephasing. Considerations of this kind are particularly 
relevant for quantum-information processing, where one 
necessarily deals with nonequilibrium situations involving 
excited qubit states of finite energy. The situation is differ- 
ent for the weak-localization problem (and similar trans- 
port interference effects): There, one is interested in the 
zero-frequency limit of the system's linear response, which 
depends on the coherence properties of arbitrarily low-lying 
excited states. It is the subtleties associated with a path- 
integral description of these situations which we want to 
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address in this work. 

The article is organized as follows: After a brief review 
of the influence functional and its meaning in semiclassi- 
cal situations (sec. II), we will rewrite the expressions for 
Si and the Golden Rule decay rate, in order to compare 
the two (sec. [II). In doingj-sa. we will closely follow the 
analysis of Cohen and ImryffiJ. Then, we will show how 
and why the effects of Sr and Sj may compensate each 
other in the derivation of a decay rate starting from the 
path integral expression for the time-evolution of the den- 
sity matrix (sec. IV), even though they cannot cancel each 
other in the influence action. There, it may be observed 
that drawing conclusions about decay directly from the ex- 
ponent iSr — Si of the influence functional is usually not 
possible. The crucial cancellation takes place at a later 
stage of the calculation, after proper integration over the 
fluctuations around the classical paths and after averaging 
over the initial state. In order to prove that this compen- 
sation takes place not only in lowest order perturbation 
theory, we will then specialize to exactly solvable linear 
dissipative systems, i.e. the damped harmonic oscillator 
(sec. [v|) and the free particle (sec. VI). We will also point 
out that there is an important difference between the os- 
cillator and the free particle: Dropping Sr has a much 
more drastic effect on the former, leading to an artificial 
finite decay rate of the ground state at zero temperature. 
However, in order to discuss the importance of Sr for the 
calculation of dephasing rates (involving decay of excited 
states ), w e have to extend our analysis to nonlinear models 
(sec. VII). We will explain in a more qualitative fashion 
why the essential insights gained from the exactly solvable 
models should remain valid both for the nonlinear models 
as well as for systems of degenerate fermions. 



the Feynman- Vernon influence functionaltra. Both Sr and 
Si are real-valued functionals that depend on both paths 
simultaneously. The influence functional is the overlap of 
bath states which have time-evolved out of the initial bath 
state xo under the action of either q > (-) or q < (-): 



JSr, 



(3) 



At zero temperature, the state xo is the ground state of 
the unperturbed bath, while at finite temperatures an addi- 
tional thermal average over xpiihas to be performed. From 
this representation, it followsaH that Sr, changes sign on 
interchanging q > and g < while Si always remains nonneg- 
ative - the magnitude of the overlap can only be decreased 
compared with its initial value of one. 

The meaning of Sr and Si becomes particularly trans- 
parent for an interference setup, where two wave packets 
, F> and , 5< travel along two different paths q^(-) and 
lci(')- We want to assume a situation which can be de- 
scribed semiclassically, i.e. the wave-length is supposed to 
be much smaller than the size of a wave packet and this 
again is much smaller than the typical separation between 
the paths and the typical dimensions over which the exter- 
nal potential changes. Then, it suffices to evaluate Sr and 
Si just for the combination of these two paths, since the 
fluctuations around them are comparatively unimportant. 
The environment will affect the interference pattern on the 
screen mainly by changing the interference termcl: 



p(x,x,t) « |<MM)| +|*<(M)r + 
2 Re [*>(a;,t)^<(a;,i)e 4Si? ^ 



(4) 



II. THE INFLUENCE FUNCTIONAL AND 
DEPHASING IN SIMPLE SITUATIONS 

We are interested in the time-evolution of the reduced 
density matrix of a system with coordinate q which is cou- 
pled to some environment (the bath). If system and bath 
had been uncoupled prior to t = 0, then the density matrix 
at a later time t is linearly related to that at time 0: 

p(<it > <it > *) = J <ko <ko J (<it > if ko,qo; t)p{qo , Qo > °) • 

(i) 

The propagator J on the right-hand-side of this equation 
is given by 

J = J J Dq>Dq< exp [i(£> - S<)] exp [iS R - Si] . (2) 

Here we have set % = 1 and introduced an abbrevi- 
ated notation: The path integral extends over all "forward" 
paths q > (-) running from the given value of q^ - to q> , like- 
wise for the "backward" paths q < (-). The value of the ac- 
tion of the uncoupled system, evaluated along q >ty< \-), is 
denoted by s£ (<) . The second exponential in Eq. (H) is 



In this equation, i $> > r < \(x, t) are assumed to represent 
the unperturbed time-evolution of the wave packets. To a 
first approximation, Srj do not enter the "classical" terms 
in the first line of Eq. (Q), since SRj[q^, q^] = 0. De- 
viations from this first approximation stem from the in- 
tegration over fluctuations away from <j£j and describe, 
for example, mass- and potential renormalization as well as 
slowing-down of the wave packet due to friction. Obviously, 
Si = Si[q^ t , qfj\ determines directly the decrease in visibil- 
ity of the interference pattern while Sr = Sr^^, qfj] only 
gives a phase-shift. Averaging the interference pattern over 
different configurations of the external potential (appropri- 
ate for impurity averaging in mesoscopic samples) may fur- 
ther decrease the visibility. However, if the two paths are 
time-reversed copies of eaah other (as is the case in discus- 
sions of weak-localizationccl) , the phase difference between 
i & > and , F< vanishes in a time-reversal invariant situation, 
so that the impurity-average over the corresponding phase 
factor does not lead to a suppression in the situation with- 
out the bath. On the other hand, the average over exp^Sfj) 
will, in general —decrease further the magnitude of the in- 
terference termtil. In any case, given the simple physical 
picture presented here, one would not expect Si and Sr to 
be able to cancel each other's effects, since they represent, 
respectively, the real and imaginary part of an exponent. 
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For the special case of a force F deriving from a bath of 
harmonic oscillators (linear bath), with vanishing average 

F ) = and a linear interaction V = —qF, we have: 



Sl = L dtl Io 1<lt2 ~ ^ RC ^ ~ ^ (5) 

Sr = - fdtJ dt 2 (q> - q<) Im (f x F 2 ) + 9 2 < )(6) 



o Jo 



Here we have used the notation q> = q > (ti). The an- 
gular brackets denote averaging over the equilibrium state 
of the unperturbed bath. Si only depends on the sym- 
metrized part of the bath correlator Re (P{t)F(0)) = 

($F(t), F(0)|\ /2, which becomes the classical correla- 
tor (F(t)F(O)} for the case of classical Gaussian random 
noise. In the latter case, Sr vanishes and the influence 
functional exp(— Si) is simply the classical average of the 
phase factor 



exp 



l , (q<(T)-q>(r))F(T)dT 




(7) 



If one drops Sr in a dephasing calculation (noting that 
it does not enter dephasing in semiclassical situations such 
as those that can be described by Eq. (Q)), one effectively 
replaces the quantum bath by a classical fluctuating force 
whose correlator is determined by the symmetrized part 
of the quantum correlator. This contains the zero-point 
fluctuations, since 



(2n(uj) + 1) cos(o;t) 



{Q(t),Q(0)} 
coth(w/2T) cos(wr 



for the coordinate Q of a single bath oscillator of frequency 
w, with n{uj) being the Bose distribution function that van- 
ishes at T = Q. 

Although replacing the quantum bath by a classical noise 
force seems to be a drastic step, it can lead to correct pre- 
dictions for the dephasing rate in semiclassical situations 
such as the one discussed above, even at zero tempera- 
ture. It is instructive to observe how this comes about 
in an exactly solvable model. A particularly simple, situ- 
ation is the one analyzed by Caldeira and Leggettll3 (see 
also Ref. [TJ]). They considered a damped quantum har- 
monic oscillator, where the initial state consisted of a su- 
perposition of Gaussian wavepackets, one centered at the 
origin, the other at a distance z. In the course of time, 
the displaced wavepacket oscillates back and forth in the 
oscillator potential well. Whenever the packets overlap, an 
interference pattern results (due to the difference in the re- 
spective momenta) . The environment leads both to damp- 
ing and dephasing, where the latter typically proceeds at 
a much faster rate. For our purposes, we are interested 
in the limit of small damping, where, to a first approxi- 
mation, the center-of-mass motion of the wave packets is 



not appreciably altered by friction in the period of the os- 
cillation. In this case, it turns out that, indeed, the re- 
sult predicted by the approximation Eq. (Q) for the at- 
tenuation of the interference pattern is correct. This can 
be seen in Ref. |l5| by taking the limit of weak damping 
(7 — > 0, uj — > u>r) in their exact result for the exponent of 
the attenuation factor (see Eq. (2.13) of Ref. |l5|) and com- 
paring this to Si (see Eq. (JEj)), which is to be evaluated for 
the pair of classical paths followed by the two wave packets, 
q^l (t) — z cos(ujt) and = 0. In terms of the quantity 

Coo listed in appendix [A[ of the present work, both results 
may be obtained by multiplying Coo by sin(J)t) 2 and setting 
7 = 0, Co = uj. In addition, in order to obtain Si[q^, q^j], 
the factors sin(a;(f — t\)) and sin(w(t — t 2 )) must be re- 
placed by cos(wii) and cos(wi 2 ). The results obviously co- 
incide for the points in time when the wave-packets meet 
(ujt = ir/2 + mr). We emphasize that the correspondence to 
the semiclassical result holds only for a situation far away 
from the ground state of the system, where at least one 
of the wave packets is in a superposition of highly excited 
oscillator states. 

At zero bath temperature, the dephasing is purely due 
to spontaneous emission, as pointed out in the discussion 
of Ref. |l5|. Energy is transferred from the oscillating wave- 
packet into the bath and the system relaxes towards lower- 
energy states. The dephasing rate is proportional to the 
total rate— r o „t at which the system leaves a given en- 
ergy levelta Of course, this rate may be obtained using 
the Golden Rule (i.e. a master equation description) only 
for weak-coupling, but the qualitative picture seems to be 
generalM In the correct description of the physical situa- 
tion considered here, the rate T out is given entirely by the 
rate of spontaneous emission, T^™. However, if Sr is ne- 
glected, then Si describes a classical noise force (equivalent 
to the zero-point fluctuations of the bath at T — 0) and 
there will be both induced emission and absorption, pro- 
(8) ceeding at equal rates = Tf^ d . Therefore, the system 
is also excited by the bath in that approximation. Never- 
theless, we have pointed out above that the total dephasing 
rate comes out right. The reason is the following: The rate 
T out is the same in both cases, because the rate of spon- 
taneous emission in the correct description is exactly twice 
that of induced emission in the approximation: 



sp 



-pem 1 j^abs 
ind ' ind ' 



(9) 



At this point, it is easy to see the physical reason why such 
an approximation may fail near the ground state of the 
system. Then, the transitions downwards in energy may 
be blocked, which completely suppresses T^™, but not Tf^ d . 
We will make this argument more precise in the following 
sections. 



III. DECAY RATES FROM Si AND GOLDEN 
RULE: DEPENDENCE ON THE SPECTRA OF 
BATH AND SYSTEM MOTION 

The growth of Si with time depends on the spectral den- 
sity of the bath fluctuations at frequencies which appear in 
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the system's motion. Formally, this can be seen by intro- 
ducing the spectrum of the system motion related to the 
given pair of paths q >< - <s> (following the ideas of Refs. [s|]l2]): 

P{u,a) = ±- [ + °° dre^{ q >{s+ T -)- q <{s+ T -)) 
2tt J_ oq 2 2 

Ad > (s- T -)-q < {s- T -)) (10) 

Here, s = {t\ + *2)/2 and r = t\ — t-i are sum and dif- 
ference times. Therefore, P is the Fourier transform with 
respect to r of the g-dependent terms in Si (Eq. (H)). We 
takegXO^') to be zero whenever t' falls outside the range 
[0,t]. 

Furthermore, we define {FF)^ to be the Fourier- 
transform of the symmetrized correlator of F, 

{FF) U ee i- J*"dTe*"({F(T),F(0)}) , (11) 

so it is real and symmetric in frequency. The same holds 
for the system spectrum P(ui,s). 

Using this, Si can be expressed in the following way: 

pt r+oo 

Si = n ds du (FF) u P(-cj,s). (12) 

J0 J -oo 

If the dependence of P on s is not essential, then Si 
grows linearly with time t, at a rate given by the "overlap 
of bath and system spectra". Similar expressions caitrbe 
derived for a spatially inhomogeneous fluctuating forceStj, 
leading to an additional fc-dependence. 

Applying this kia-d of reasoning to a single electron mov- 
ing in a dirty metalB, the system motion is found to contain 
frequencies up to (at least) 1 /r e ; , such that the growth of Si 
with time depends on the bath-spectrum up to this cutoff, 
including the zero-point fluctuations of the corresponding 
high-frequency bath modes (which become important at 
low temperatures). As explained above (sec. II), the de- 
scription of dephasing in terms of Si alone may be trusted 
in a semiclassical situation, where the electron is in a highly 
excited state. This holds even at zero temperature, when 
the qualitative physical picture is essentially the same as 
in the model of oscillating wave packets due to Caldeira 
and Leggett, discussed above: Dephasing is due to sponta- 
neous emission of energy into the bath. The contribution of 
frequencies up to l/r e ; then implies that the spontaneous 
emission (and the resulting dephasing) is facilitated by the 
impurity scattering. The physics behind this is well-known 
in another context: In quantum electrodynamics, the emit- 
ted radiation would be called "bremsstrahlung", since it is 
the scattering off an external potential that induces the 
electron to emit radiation. 

On the other hand, one can express in a similar manner 
decay rates from a simple Golden Rule calculation^. For 
the decay of an initial state \i) , we have 

r, = 2ir I dto (pF) u (qq) { ll ■ (13) 



Here \FF^ and (<?<?) ^ are the Fourier transforms of 

the unsymmetrized correlators of F and q, taken in the 
equilibrium state of the bath or the initial state of the 
system, respectively: 

(FF) u = ±f*re*"(fiiT)fi(0)) 

= dre^(i\q(r)q(0)\i) (14) 

The important point to notice is that, near the ground 
state of the system, the quantum correlator (<?<?) |y is very 
asymmetric in frequency space, as the system can mostly 
be excited only (ui > 0). At low temperatures, the same 
holds for the bath. Thus, the decay rate fll3|), contain- 
ing the product of correlators evaluated at u> and — cj, is 
very much suppressed below the value that it would acquire 
if either the bath-correlator or the system-correlator were 
symmetrized (becoming symmetric both in time and fre- 
quency). Since dropping Sr is equivalent to symmetrizing 
the bath correlator and, furthermore, semiclassical calcula- 
tions give a symmetric spectrum P(u>, s) of the system mo- 
tion as well, it becomes clear why there are situations when 
the decay rate, as deduced from Si, is finite at low tem- 
peratures, while the Golden-rule decay rate vanishes. The 
question then arises whether any procedure that amounts 
to symmetrization of the correlators, thus leading to dras- 
tically wrong results for Golden Rule decay rates at zero 
temperature, may be justified to discuss dephasing using 
a path-integral approach. Observations such as this have 
led Cohen and Imry to conclude that, in their semiclas- 
sical analysis of electron dephasing inside a meta£3, the 
contribution of the bath's zero-point fluctuations to the 
dephasing rate should be dropped. Their argument was 
not a mathematical proof, but rather drawn from physi- 
cal intuition and analogies. In the following sections, we 
try to elucidate the importance of Sr in descriptions of 
dephasing and decay near the ground state of a system, 
demonstrating exactly how a cancellation between Sr and 
Si may arise. 



IV. CANCELLATION OF S R AND Si IN 
LOWEST-ORDER PERTURBATION THEORY 

We consider a system which is in its ground state at t = 
before being coupled to a linear bath. From the Golden 
Rule, we expect there to be no finite decay-rate for the 
ground state at zero-temperature. It is well-knowna how 
to derive a master equation from the influence functional, 
by expanding it to lowest order in the exponent, iSu — Si. 
We display that calculation here in order to point out where 
and how Sr and Si do cancel. The time-evolution of the 
probability to find the system in its unperturbed ground 
state is: 



P {t) = (tf |#i)|tto> 
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= / dq>dqf Po {qf,q>)p{q>,q< 1 t) 

= J Dq>Dq< PQ ( q <,q>)e^- S o^ S «- S ' Po ( q >,q<) 

« 1 + J Dq>Dq<e i ^- s o)p (q<, q >) 

x(iS R -Si)p (q>,q<). (15) 

Here po is the density matrix of the unperturbed ground 
state, and the path integration also includes integrals over 
the initial and final coordinates. Now one can insert the 
expressions for Sr and Si from Eq. @. The integration 
over the trajectories and the endpoints produces the (un- 

I 



symmetrized) correlators of the system coordinate, such as 
(for ti > t 2 ): 



Dq> e iS o %(q>)q>( tl )q> (i 2 )tt„ (q> ) 



(*o|g(*i)?(*2)|*o)e- 



-iEnt 



(16) 



Evidently, both the unperturbed action Sq and the inte- 
grations involving the state are essential for obtaining 
the correct correlator. In the long-time limit, we obtain the 
decay rate which also follows from the Golden Rule (|l3|) : 



ft 

P (t) ml- I dh dt 2 2Re 



o 







F(ti)f(t a ))(*o|$(*i)?(ta)|*o> 



1 - Tot . 



(17) 



Only the term growing linearly in time has been retained 
at the end of this equation. 

Note that, of course, there will always be a small reduc- 
tion in the probability of finding the unperturbed ground 
state after switching on the interaction, since the ground 
state of the coupled system contains contributions from 
other system states as well. This point has been discussed 
in more detail recently for the case of the damped harmonic 
OSCillatorHL It does not contribute to r which describes 
the decay linear in time. 

We emphasize that the decay of Po (*) is not governed by 
Si alone: If Sr were dropped, then only the (real-valued) 
symmetrized version of the bath correlator would appear 
m Eq. (0). In that case, the decay rate would only depend 
on the (real- valued) symmetric part of the system correla- 
tor as well. Therefore, a finite decay rate of the ground 
state would result even at zero temperature, where really 
there would have been none. The physical reason is the 
heating introduced by the classical noise field, whose cor- 
relator contains the zero-point fluctuations of the original 
quantum bath, as we have discussed before. 

On a formal level, we may argue that the decay of 
the density matrix is not governed by Si alone since the 
"weighting factor" that is used when "averaging" over many 
paths contains the phase factor exp(i(S' > — Sq)). There- 
fore, such an average is not the same as a classical aver- 
age, for which the decay could not be overestimated by 
dropping Sr, since then we could use \ (exp(iSR — Si))\ < 
(exp (—Si)). Furthermore, the integration over the density 
matrix of the initial state is obviously essential, as it is 
needed to produce the correct form of the system correla- 
tor. Both facts mean that it would be premature to draw 
any conclusions about dephasing and decay at an early 
stage of the calculation, by merely looking at the influence 
action. 



V. EXACTLY SOLVABLE LINEAR SYSTEMS: 
"CANCELLATION TO ALL ORDERS" FOR THE 
DAMPED OSCILLATOR 

In order to show how Sr and Si can cancel also be- 
yond lowest-order perturbation theory, we will now turn 
to linear quantum dissipative systems. Although the ex- 
act solutions for the damped harmonic oscillator and the 
free particle have been well-known for a long timetj, we 
will review the essential steps in the derivation, pointing 
out where Sr and Si do enter. We will first turn to the 
oscillator, for which the Golden Rule result suggests that 
one would obtain an artificial decay of the ground state at 
zero temperature if Sr were neglected. For simplicity, our 
discussion is restricted to T = 0, since that is the limit 
where these effects show up most clearly. 

For any linear system which is linearly coupled to a 
bath of oscillators, the propagator J for the density matrix 
(Eq. (||)) can be evaluated exactly, since the integration 
over system paths is Gaussian. In fact, J is found to be 
given by the "semiclassical" result, i.e. an exponential con- 
taining the action along stationary paths, multiplied by a 
prefactor which does not contain the endpoints of the paths 
(a specialty of linear systems): 



•Agf^ko ,«o ;*) 



1 



N(t) 



JS[q> v q<] 



(18) 



with S = Sq + Sq + Sr + iSi. The paths <j£j make 
the full action stationary, 



SS 



^ = 

< u ' 



(19) 



and fulfill boundary conditions of the form q^ (0) 
The prefactor N(t) can be obtained most easily from the 
condition for normalization of the density matrix, 
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dq t J(qt,qt\qo,Qo;t) = 1 . 



(20) 



N(t) is found to be independent of Si, but it does depend 
on Sr (see, for example, the general proof in app. E of 
Ref.f). 

In the following, we will turn to the special case of the 
Ohmic bath, which leads to a velocity-proportional fric- 
tion force and has a power spectrum rising linearly at low 
frequencies (for T = 0) : 



2tt 



(21) 



We will argue below that, in the case of the damped 
harmonic oscillator, no essential qualitative result will be 
changed by using other bath spectra, unless these have an 
excitation gap exceeding the oscillator frequency. 

As usual, we introduce the center-of-mass coordinate 
(q > ( T ) + < Z < ( T ))/2 and the difference coordinate 



R(r) 

r(r) ee 9 >(r)-£( 
obtained from 



r) in order to write down the equations 
for the case of the Ohmic bath: 



dr 2 



dr 2 



= 



(22) 



d 2 R dR 2 „ 



i\ dr' Re{F(T)F(T')jr(T'X23) 



Here 7 ee rj/m is the damping rate and ujq is the 
unperturbed frequency of the oscillator. Note that the 
second equation leads to a complex-valued solution R(-). 
It is also p-Bpasible to formulate the calculation slightly 
differentlytSEj, by using stationary solutions with respect 
to the real part of S only. In any case, inserting the so- 
lutions R (= R c i) and r (= r c i) into the action S shows 
that the imaginary part of the action S c i = is 
determined directly only by Si. As expected, S c i turns out 
to be a bilinear expression in the endpoints R t , Rq, r t , r : 



Re Sri = 



ImSrj, = 



Sr — 

RtnLtt + Rtr L m + R r t L ot + i? r £oo(24) 
Si = r 2 C m + 2r Q r t C ot + r 2 t C tt . (25) 



All the coefficients are real- valued functions of time. For 
our purposes, it is sufficient to know that the entries of 
the matrix L are independent of Si (while they do depend 
on the damping 7 from Sr). In contrast, the entries of C 
depend both on Si and Sr. The dependence on Sr arises 
only because the path r, which is inserted into Eq. (|5), 
is affected by the friction described by Sr, see Eq. (|22). 
The entries of C are proportional to the strength of the 
bath fluctuations (i.e. the magnitude aL(Fi 7 ') w ). Explicit 
expressions for all of these coefficients 2 ^ can be found in 
Ref. [l8| (cf. their section 6). For reference, the quanti- 
ties used in the following discussion have been listed in 
appendix [A|. 



At this point, evaluation of the coefficients for ImS c i 
shows that they will grow in time beyond all bounds, re- 
gardless of whether the path r is calculated by takin g in to 
account Sr or neglecting it (setting 7 = in Eq. (22)). 
However, no conclusion about dephasing and decay can be 
drawn from this, since the foregoing discussions lead us to 
expect that any potential cancellation between Sr and Si 
will take place only after proper integration over the initial 
density matrix. This can be deduced from the derivation of 
the Golden Rule decay rate given in the preceding section. 
It is also very reasonable that Si grows without bounds for 
any semiclassical path r(-), since all paths contribute more 
or less to the time-evolution of all eigenstates, and decay 
of excited states is perfectly correct in the model of the 
damped harmonic oscillator. 

The initial ground state density matrix has the form 



Po(<?o- 



«o < ) 



cxexp(-( (?0 > 2 + (Zo < 2 )/(4( ( z 2 ) )) 
r 2 /4)/(2(«f) )) , 



exp 



(26) 



with (q 2 ) = l/(2muo) as the square of the ground state 
width. After performing the Gaussian integrals over Rq 
and r , which are needed to evaluate Eq. (|j), we obtain 
the final result for the density matrix p(q> , qf , t) at a time 
t after switching on the interaction with the bath. It has 
the general form 



N(t) 



exp (-(aR 2 + bR t r t + cr 2 )) 



(27) 



In particular, the prefactor of R 2 gives the width of the 



probability distribution at time t, 



(q 2 (t)) = (2a) 



(28) 



Let us now look at the behaviour of the width in order to 
analyze the effect of dropping Sr in the calculation. The 
expression iaL^(q 2 (t)) has been derived already by Caldeira 
and LeggettES (cf. their Eq. 6.34): 



(q 2 (t)) = (2C7 00 + muj /2 + LlJ{2moj Q )) /L 2 t0 . (29) 

This expression clearly contains quantities both from the 
imaginary part (Coo) and from the real part (Lqo, L t o) of 
the action. Formally speaking, real and imaginary parts 
have become intermixed when integrating over Rq, due to 
quadratic completion of the expression zi?o(r t Lo t +roLqo) — 
mwoi?o in the exponent (the former term is from Eq. (|24|), 
while the latter term is from the initial density matrix, Eq. 
(p6|)). Physically, this result is to be expected, since the 
growth of the width of the probability distribution will be 
governed by the balance between the bath fluctuations (Si) 
and the friction (Sr). 

The complete time-evolution of (q 2 (t)) after switching 
on the interaction at t — is displayed in Fig. [j]. This 
includes cases where the value of the damping constant 7 



has been artificially set to zero or to other values different 
from that prescribed by the fluctuation-dissipation theorem 
(FDT), which connects 7 to the strength of the bath fluc- 
tuations described by (FF) U (see Eq. ( |2l| ) for the case of 
T = 0). If Sr is neglected completely (no friction, 7 = 0), 
heating takes place, causing the variance to grow linearly in 
time, at a rate given by the Golden Rule expression involv- 
ing the symmetrized bath correlator. Formally, the growth 
of Coo with time cannot be compensated in that case, since 
Loo and L t o acquire their original unperturbed values. For 
small 7 > 0, the width saturates at a value much larger 
than that of the ground state. On the other hand, if the 
friction is too strong, the width saturates at a value below 
that of the ground state, and the Heisenberg uncertainty 
relation is violated (since (p 2 ) alsojshrinks) . This point 
is discussed in the book of MilonniEil in the context of an 
atom interacting with the vacuum electromagnetic field. If, 
however, 7 has the correct value prescribed by the FDT, 
(q 2 ) changes only slightly from its unperturbed value. 

In fact, for a linear system like the quantum harmonic 
oscillator, its behaviour under the action of the linear bath 
can be [described and understood entirely using a classical 
pictureO. The reason for this is as follows: one starts out 
with an initial state for system and bath that has a pos- 
itive definite Gaussian Wigner-density, which can be in- 
terpreted directly as a classical phase space density whose 
time-evolution then corresponds one-to-one to the evolu- 
tion of the Wigner-density (see appendix [b]). Thus, the 
quantum dissipative dynamics may be described by a clas- 
sical Langevin equation including friction and a fluctuating 
driving force: 

d 2 q dq 9 1 , „ , „ 

In this equation, the damping term incorporates the ef- 
fects of Sr. It counteracts the fluctuations of F(-) that 
are described by Si in the path-integral picture and whose 
correlation function is, therefore, given by the symmetrized 
part of the bath correlator. 

The behaviour of the width has been discussed above 
using the influence functional, but it can be understood 
more easily in this picture. We can use the susceptibility 
of the damped oscillator to find its response to the force 
fluctuations F, which, at zero temperature, are entirely due 
to the zero-point fluctuations. In particular, for the limit 
t — » 00, we obtain 




Q| , 1 , 1 , 1 , 1 

10 20 30 

t 

Figure 1: The variance a 2 (t) = (q 2 {t)) of the probability 
density p(q,q,t) for the damped harmonic oscillator, plotted 
in units of the unperturbed ground state variance a 2 , as a 
function of time t for the following friction strengths (top to 
bottom): 7/(27r (FF) w=1 ) = 0, 0.5, 1, 2 (other parameters are 
m = Uq = 1). A ratio of 1 corresponds to the correct value pre- 
scribed by the FDT, while 7 = means Sr has been neglected 
completely, such that cr 2 {t) grows without bound. The dashed 
lines give the limit for t — > 00 (see Eq. dil])). 



correct equilibrium state (obtained by fully keeping Sr and 
Si), the width at t — > 00 can also be found by applying the 
FDT to the oscillator coordinate (cf. Eq. 6.37 of Ref. ||). 
In contrast, the effect of dropping Sr, is obtained by letting 
7 -> in Eq. ©, while {FF) U is kept fixed. 

The example of the damped harmonic oscillator demon- 
strates that the proper behaviour near the ground state of 
the system cannot be observed at the early stages of the 
calculation, by looking at the action evaluated along clas- 
sical paths (Eqs. (pi|), (p5|)), regardless of whether these 
paths properly include the damping or not. Only after 
correct integration over the initial state density matrix the 
contributions from Sr and Si can compensate each other 
(in Eq. (|29|)), such that no artificial finite decay rate of the 
ground state results, if Sr is kept. 



+00 



where 



duj |xH| 2 (FF) U 



1/m 

- UJ 2 — 1077 



(31) 



(32) 



is the susceptibility of the damped oscillator, and (FF)^ is 
the symmetrized power spectrum (see Eq. (pi])). For the 



VI. POWER-LAW BEHAVIOUR IN QUANTUM 
BROWNIAN MOTION 



We now turn to a discussion of the free particle, i.e. 
the Caldeira-Leggett model of quantum Brownian motion. 
There, the reasoning is the same as before in principle: 
Neglecting Sr means that the zero-point fluctuations of the 
bath contained in Si may heat up the particle. However, 
there is an important difference: The spectrum of the free 
motion does not contain any resonance peak, unlike that of 
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the harmonic oscillator. Therefore, only the low-frequency 
components (uj — » 0) of the bath spectrum will contribute 
to heating and decay (compare E«h|-(t2l|), which results 
in a peculiar long-time behaviounisllialll that cannot be 
captured using the Golden Rule, as has been emphasized 
in Ref. [| 

One can observe the distinction between oscillator and 
free particle most easily in the picture of the classical 
Langevin equation ( |3C| ) introduced above. In appendix ||, 
it is shown how the density matrix propagator J may be ob- 
tained in general via the time-evolution of the Wigner den- 
sity. This approach is physically more transparent than the 
path-integral calculation. For the purposes of our present 
discussion, however, it will suffice to consider the time- 
evolution of the diagonal elements of the density matrix in 
momentum space. These have also been analyzed in Ref. 0, 
in order to demonstrate the failure of the Golden Rule cal- 
culation to describe the decay of the original ground state 
at zero temperature. We will confirm and explain the out- 
come of this analysis using a simple argument based on the 
Langevin equation (30). Then we will point out the role of 
Sr and describe the important difference to the damped 
oscillator. 

If the particle has momentum po at time 0, its momen- 
tum p at time t will be determined by the fluctuating force 
F(-) in the following way (by solving Eq. (|30|) for ujq = 0): 



p = p e-" /t + / dse-^^Fis). 
Jo 



(33) 



Since F(-) is a Gaussian random process, the probability 
density of finding a momentum p at time t is a Gaussian, 
of variance 



(Sp 2 ) = / dsi 
Jo 



ds 2 e-^ 2t - s ^ (F( Sl )F(s 2 )) . (34) 



The behaviour of (Sp 2 ) as a function of time is as follows 
(at T = 0): As long as t -C I/7, friction is unimportant and 
it is the double time-integral over the correlator of F which 
is to be evaluated. For the Ohmic bath considered here, 
the power spectrum is relatively strong at low frequencies 
(rising linearly with u, see Eq. (^l|)). This leads to a 
slow decay in time, (F(t)F(0)) oc l/t 2 , which results in a 
logarithmic growth of (Sp 2 ), 



(5p 2 ) 



2-ln(w c t) 

7T 



(35) 



Here uj c is the cutoff frequency of the bath spectrum and 
we have assumed u> c t ^> 1. It is this logarithmic behaviour 
which cannot be obtained using the Golden Rule approx- 
imation, to be discussed below. At later times, t ;» 1/7, 
the growth saturates at a constant value, 



(<5p 2 >^ln(a, c / 7 ). 



(36) 



The perturbation expansion in the coupling between sys- 
tem and bath can be carried out by expanding the influence 



action in the propagator J (Eq. (|2|)) in powers of iSr — Si 
see Eq. ( fl5| ) in the present article or Eq. (E13) of Ref. 

We want to discuss the time-evolution of the diagonal 
elements p pp (t) of the density matrix in momentum space. 
The evolution starts from the equilibrium density matrix, 
which is a Gaussian momentum distribution of variance 
(po) = f>T (the following formulas will be analyzed for 
arbitrary T). After a time t, the density matrix is still a 
Gaussian, but of variance (compare Eq. (P3[)): 



(P 2 ) = (PI)< 



(Sp 2 ) 



(37) 



From this we can easily derive the result for Sp p p\t) 
which could alternatively be obtained by first expanding 
J to first order in %Sr — Si, going over to the momentum 
representation and finally integrating over the initial den- 
sity matrix (as was done in Ref. 0; see Eq. (E18), which 
is, however, still written in terms of matrix elements of the 
coordinate operator): 



p PP {t) 



exp 



2(p3) 



Mo) • (i 



(p 2 ) (pI) ( P 



2<P§> \{PI) 



(38) 



In evaluating (p 2 ), only the terms up to first order in 
7 and (FF)^ must be kept for the purposes of this ex- 
pansion. At finite temperature T, the behaviour of (p 2 ) 
at long times t 3> 1/T is governed by the linear decrease 
stemming from (p 2 ,) e~ 2lt re (p 2 ,) (1 — 2 7 t) and the linear 
increase from (Sp 2 ) ~ 2-qTt. If the initial density matrix 
really describes the equilibrium distribution, then detailed 
balance holds and these terms cancel. For times (<1/T 
(or arbitrary times at T — 0), the time-evolution of (p 2 ) 
is governed by (Sp 2 ), evaluated for the correlator {FF) U 
taken at T = (which only contains the zero-point fluctua- 
tions of the bath). For consistency of the expansion, (Sp 2 ) 
in Eq. (Q) has to be evaluated by setting 7 = 0. There- 
fore, only the discussion given above for times t <C 1/7 
turns out to be relevant. The logarithmic growth of (Sp 2 ) 
given in Eq. (|35| ) corresponds to a power-law behaviour of 
the (exact) time-evolution of p pp (t) for t <C I/7. This log- 
arithm, that appears in the perturbation expansion (|3q), 
is "overlooked" in the Golden Rule approximation, where 
only terms growing linearly with time are .kept. Therefore, 
the Golden Rule rate turns out to vanishEI. These power- 
laws are characteristic of the density matrix propagator 
J of quantum Brownian motion at zero temperature (see 
appendix [§|) . 

On the other hand, the effect of neglecting Sr may be 
analyzed by setting the friction constant 7 to zero in Eq. 
(|5]). Then, the spread (Sp 2 ) of the momentum distribu- 
tion grows without bounds, which is obviously not the cor- 
rect physical behaviour. It is qualitatively similar to the 
heating produced in the damped oscillator model. How- 
ever, since the growth proceeds only logarithmically in 
time, the artefacts of this approximation, when applied 
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to the free particle, are not nearly as drastic as the finite 
decay rate of the ground state observed for the oscillator. 
In particular, within perturbation theory, no qualitative 
change is obtained by dropping Sr for the free particle (no 
finite decay rate is produced). 

Furthermore, if we had considered a super-Ohmic bath 
(e.g. due to phonons), whose spectrum decays faster than 
uj 1 for u> — * 0, the growth of (p 2 ) would saturate even 
without friction (as it should, since there is no velocity- 
dependent friction for such a bath). In contrast, the be- 
haviour of the harmonic oscillator discussed above would 
remain qualitatively the same as for the Ohmic bath, since 
its decay depends on the bath spectrum at the resonance 
frequency of the oscillator and is not affected in any essen- 
tial way by the details of the spectrum at low frequencies. 
This is consistent with the fact that in the oscillator case 
already the Golden Rule is sufficient to obtain an essen- 
tially correct picture of the artificial decay produced by 
dropping Sr, while, for quantum Brownian motion, it fails 
to describe the subtle power-laws associated with the low- 
frequency (long-time) properties of the Ohmic bath at zero 
temperature. 

The difference between oscillator and free particle also 
shows up clearly in the behaviour of the imaginary part of 
the influence action: For the free particle, Si always grows 
only logarithmically with time t (at T = 0), independent 
of whether one keeps Sr or sets it to zero (i.e. 7 = in 
the equations of motion). This can be derived from the 
formulas given in appendix [b] or those of Ref. 0, Appendix 
E (or Ref. |l8[ section 6, with the proper limit Wfl->0 for 
the free particle). 

Hence, we conclude that for the oscillator (or any motion 
containing an extended spectrum) the effects of dropping 
Sr are much more pronounced than for the free motion, 
since they lead to an artificial finite decay rate. We have 
also observed that this effect can already be understood 
within the framework of the Golden Rule approximation, 
whereas that approximation is incapable of describing the 
more subtle "power-law decay" found for the free particle. 



VII. QUALITATIVE DISCUSSION OF OTHER 
MODELS: NONLINEAR COUPLING AND PAULI 
PRINCIPLE 

It could be argued that our discussion of the damped 
oscillator has only demonstrated the importance of Sr in 
preventing an artificial decay of the ground state and has 
little to do with dephasing. Indeed, if one wants to discuss 
dephasing, one should rather consider the decay of a co- 
herent superposition of the ground state and some excited 
state. This is relevant both for arbitrary nonequilibrium 
situations as well as for the calculation of the system's lin- 
ear response, where the perturbation creates a superposi- 
tion of excited states and the ground state. In the example 
of the damped harmonic oscillator, full relaxation into the 
ground state will always take place (at T = 0) , for any such 
excited state. This holds regardless of the bath spectrum, 
provided the latter does not vanish around the resonance 
frequency of the oscillator. For weak coupling, the decay 



of the density matrix (and, therefore, the linear response) 
may be described conveniently using a simple master equa- 
tion approach. What is more, even if one neglects Sr in 
a path-integral calculation, the predicted decay rates will 
be correct, excepting only that for the ground state, which 
we have discussed above. The reason for this, however, 
is a specialty of the linearly damped harmonic oscillator. 
Since the bath couples to the coordinate operator, whose 
matrix elements only connect adjacent oscillator levels, the 
relevant "system correlator" (in Eq. jlO|)) turns out to be 
symmetric for any state but the lowest one. Therefore, the 
decay rate will depend only on the symmetric part of the 
bath correlator as well, see Eq. JTij). Hence, for all the 
excited states, the evaluation of the total decay rate is not 
affected by dropping Sr. This point has already been dis- 
cussed in connection with the model of two oscillating wave 
packets, at the end of section II. 

Since we are interested in demonstrating the importance 
of Sr for the dephasing of a superposition of low-lying lev- 
els, the natural choice is to consider a bath with an exci- 
tation gap larger than the resonance frequency of the os- 
cillator (see Fig. ||). In that case, we expect the lowest 
levels to remain coherent in the correct description. Un- 
fortunately, in the model of the linearly damped oscillator 
discussed above, we would not obtain any finite decay rate, 
regardless of whether Sr is kept or not. This is because 
the corresponding classical noise force is out of resonance 
and cannot heat the oscillator. 

Therefore, in order to analyze a situation where the 
importance of retaining Sr is displayed even for excited 
states, we have to go beyond exactly solvable linear mod- 
els. As a consequence, the following discussion is neces- 
sarily incomplete insofar as we cannot give exact proofs of 
the statements to be made below. These statements will 
be based on the experience acquired in the simpler mod- 
els discussed above. We will also make use of some Golden 
Rule type arguments, noting that the Golden Rule has been 
sufficient to understand the behaviour of the damped oscil- 
lator in a qualitatively correct way. In any case, we believe 
it is useful to contrast the results obtained by dropping Sr 
with the "commonly accepted" picture. 

Let us, therefore, consider the following model: We re- 
tain the harmonic oscillator, but change the coupling be- 
tween system and bath from V — —qF to V = —f(q)F, 
with a nonlinear function /. This is called state-dependent 
frictionEi In that case, the transition matrix elements of 
the system operator /(g) will, in general, be nonzero be- 
tween any two states, which is the important difference 
to the linear case. If we introduce a bath spectrum with 
a gap larger than nw , then, for sufficiently weak coupling 
and according to a simple Golden Rule calculation, the first 
n excited states will not decay at zero temperature. Any 
coherent superposition of these states will therefore remain 
coherent. The full, non-perturbative picture will be slightly 
more complicated, but the essential point should remain 
the same: The first excited states will acquire an admix- 
ture of other levels and become entangled with the states 
of the bath. They will experience some frequency shifts 
and the transition matrix elements will be renormalized. 
Of course, the new eigenstates of the coupled system will 
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Figure 2: The models described in the text. Left: A bath 
with a gap would lead to the upward transitions indicated by 
grey arrows if (and only if) Sr were neglected. Right: A sea of 
fermions may block downward transitions (dashed) even though 
they are possible in the single-particle picture. 



remain in a coherent superposition forever (by definition). 
The n+1 lowest ones still have a discrete spectrum and are 
in direct correspondence to the initial unperturbed states. 
When switching on the interaction at t — 0, a partial decay 
will result. Physically, this corresponds to the relaxation 
of the initial state into the selfconsistent coupled state of 
system and bath (see Ref . |l^ for a discussion of these issues 
in the case of the free particle). However, the initial decay 
will saturate in time on a short time-scale. During this 
transient adjustment, the reduced system density matrix 
becomes mixed to a small extent (if the coupling is weak) , 
but no non-saturated long-term decay (indicated by a fi- 
nite decay rate) results. In any case, the transient decay is 
an artefact of the procedure of suddenly switching on the 
interaction. A full calculation of, for example, the linear re- 
sponse properties of the dissipative system would start out 
with the selfconsistent ground state of the fully interacting 
system (as was done in Ref. for the free particle) . 

In contrast, dropping Sr would lead to a completely dif- 
ferent picture, since then Si would correspond to a classical 
noise force. Starting from any of the first n excited states, 
the fluctuating field would always be able to induce transi- 
tions up in energy, thus leading to a finite decay rate. The 
role of Sr in the correct description of a quantum noise 
force is precisely to cancel these upward transitions, as has 
been demonstrated in the case of the linearly damped os- 
cillator. Due to this finite decay rate, any initially coherent 
superposition would get destroyed, leading to complete de- 
phasing. No saturation of the decay could be observed, in 
spite of the gap in the bath spectrum. 

It is this behaviour that is described qualitatively cor- 
rectly by evaluating Si along semiclassical paths of the 
oscillator. The expression for Si, Eq. (||), is now to be 
changed simply by replacing q by f(q). While an (un- 
perturbed) oscillator path of the form q(r) oc cos(wo^) 
only contains frequencies ±w which cannot couple to the 
gapped bath, the function f(q{r)) will, in general, contain 
all higher harmonics as well. This directly corresponds to 
the character of the symmetrized system spectrum, i.e. the 
Fourier transform of ({f(q(r)), f(q(0))}). These frequen- 
cies contained in the system motion then will couple to the 
bath spectrum, leading to an unbounded growth of Si with 
time. 



Similarly, if we were to use a bath without a gap, but 
where the spectrum falls off fast towards low frequen- 
cies, then the correct description would also yield decay 
rates that quickly become smaller when going towards the 
ground state, while they would saturate at a finite, com- 
paratively large level if Sr were dropped. 

The present work has been motivated in part by 
the ongoing controversial debate on low-temperature 
dephasi«e-||— ,111, p- iptr— i interacting disordered fermion 
systemllHiMBBEil, which has led to a revival of in- 
terest fat, the g eneral question of dephasing in mesoscopic 
systemsM'E3t§EZfE3o. Since the usual Fey nman- Vernon 
influence functional deals only with a single-particle 
situation, it cannot be applied directly to this problem. 
However, the authors of Refs. HJt],|| have succeeded in 
deriving an extension of the usual influence functional to 
the many-fermion situation, using an exact procedure. In 
their result, the Fermi distribution (and, therefore, the 
Pauli principle) only enters Sr, while Si is unaffected 
(Eqs. (54) and (55) of Ref. ||). Therefore, the dephasing 
rate, as read off from Si, is found to equal the rate 
which would also be obtained in a purely single-particle 
calculation. Although we do not analyze the full problem 
of weak-localization here, we will use the insights gained 
above in order to explain why, in our opinion, the as- 
sumption that the dephasing rate can be derived from Si 
alone must be proven instead of being taken for granted. 
This holds even when Sr vanishes on the relevant pairs 
of classical time-reversed trajectories, since, in general, 
the integration over the fluctuations away from these 
trajectories will be essential for obtaining a cancellation 
between the effects of Sr and Si (see sec. [V), in addition 
to properly taking into account the initial density matrix. 
Furthermore, we note that in these calculations the 
detailed form of the bath spectrum at low frequencies 
turns out to be unimportant for the essential result of a 
finite zero-temperature dephasing rate. This is in marked 
contrast to the case of the-£ree particle (where the Ohmic 
bath plays a special roleca; see also Refs. p§| , f29| ), but 
similar to what is observed for the damped harmonic 
oscillator when Sr is neglected. It is for this reason that 
we have chosen a bath with an excitation gap in order to 
demonstrate the importance of Sr, since there the effects 
come out most clearly, although they are also present for 
other bath spectra (compare the discussion above). 

The model situation of the damped oscillator described 
above already contains one key ingredient related to the de- 
scription of dephasing for an electron moving inside a dis- 
ordered metal, namely an extended system spectrum: Due 
to the impurity scattering, the system spectrum (e.g., of 
the velocity operator) contains frequencies up to (at least) 
the elastic scattering rate, which may couple to the zero- 
point fluctuations contained in the bath spectrum. This 
is in contrast to the free particle, which can only couple 
to the low-frequency bath modes (see the discussion at the 
end of section |Vl| ). As before, the importance of Sr will be 
visible most easily for a bath containing an excitation gap. 
For such a bath, the free motion (without impurity scat- 
tering) will not show any nontransient decay, regardless of 
whether Sr is taken into account or not. In contrast, the 
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model to be discussed below, which is more similar to the 
situation in the disordered metal, will show a finite decay 
rate of low-lying levels if and only if Sr is neglected. 

Apart from the extended system spectrum, we have 
to take intarnStficpint another important feature of the 
calculationsa&BKL3 concerning electrons in a disordered 
metal: the use of the semiclassical analysis. For linear sys- 
tems the semiclassical result for the path- integral is exact, 
such that it can even be used near the ground state of the 
system, as we have done it here. However, in such a case, 
the size of the fluctuations around the path is comparable 
to the amplitude of the path itself, which is not the situ- 
ation in which semiclassics is usually applicable. In order 
to gain intuition for a typical semiclassical situation, we 
could reconsider the example with the nonlinear coupling 
to a gapped bath given above: Imagine an initial superposi- 
tion of fast wavepackets, whose wavelength is much smaller 
than the packet size but whose excitation energies are still 
below the bath threshold. They would show a finite decay 
or dephasing rate (for the reasons given above) , if only Si 
were used for the analysis, but not in the full calculation. 
Still, for this scenario one might argue that it has been 
clear from the outset that the semiclassical analysis can- 
not be trusted, since those high-frequency components in 
the system spectrum that are responsible for the decay are 
necessarily larger than the energetic distance to the ground 
state. 

It is only in a system of degenerate fermions that the 
following two conditions can be fulfilled all at once: On 
the one hand, the semiclassical analysis of a single (non- 
interacting) electron moving at the Fermi level is valid for 
an external potential which is sufficiently slowly varying, 
such that the system spectrum (concerning the motion of 
the single electron) only contains frequencies much smaller 
than the Fermi energy. On the other hand, the whole 
many-particle system may be near (or in) its ground state. 

While the first feature would lead one to believe that 
Sr may be omitted, we have learned from the examples 
discussed before that this is likely to be incorrect when- 
ever the second condition holds as well. We stress once 
again that the two conditions are mutually exclusive in a 
single-particle problem, which is the reason why such con- 
siderations have not played any role in influence-functional 
calculations up to now. 

In order to render the discussion concrete, we can once 
again make use of the example with a nonlinear coupling 
f(q) given above, provided we suppose the N lowest oscilla- 
tor states to be filled up with fermions initially (see Fig. ^). 
For a relatively smooth function /, the matrix elements of 
f(q) only connect states within a range much smaller than 
ep = Nu>o, which, in our example, should still be larger 
than the gap nu>o of the bath. Dephasing of a single parti- 
cle near e_p can then certainly be described fully within the 
semiclassical analysis, and, to a good approximation, us- 
ing Si alone, as has been discussed above. The same holds 
for the many-particle system, if one explicitly considers a 
classical noise force, where Sr is absent. In that case, the 
many-particle problem can really be treated as a collection 
of independent single-particle problems, as is the case for 
any external time-dependent potential. Only in the end 



an average of the full Slater determinant over all possible 
realizations of the external noise has to be carried out. 



If the dephasiH|g|rate is calculated solely from Si, it turns 
out to be finiteotdfl and not to depend at all on the distance 
to the Fermi surface. This is consistent with the fact that 
the the value of ep does not even appear in the single- 
particle calculation. On the other hand, for a quantum 
bath Sr does not vanish and should be included in the 
influence action. It is true that for the case of a highly 
excited single particle the dephasing rate comes out cor- 
rect, regardless of whether Sr is kept or not, as we have 
discussed before. However, we have also pointed out that 
there is an important physical difference between the two 
calculations. Since in the correct approach (including Sr), 
the transitions induced by the bath are purely downwards 
in energy (accompanied by spontaneous emission into the 
bath), it is reasonable to expect that, in the many-particle 
problem, they will be blocked by the Pauli principle. In 
any case, this is what is found using the Golden Rule. 
Judging from these arguments, the contribution to Si from 
the nonzero overlap of the (symmetrized) system spectrum 
with the high frequencies of the bath spectrum does not 
imply dephasing and decay, but rather a renormalization, 
similar to that obtained for an electron interacting with 
optical phonon modes (a gapped bath), leading to the for- 
mation of a polaron. The formation of the polaron will be 
visible as an initial transient decay of the single-electron 
density matrix (for the artificial case of factorized initial 
conditions), which saturates on a short time-scale. There- 
fore, to lowest order, for a bath with a gap no finite re- 
laxation rate of the lowest-lying single-particle excitations 
above the Fermi sea is expected, just as for the lower levels 
of the harmonic oscillator in the single-particle model with 
nonlinear coupling to a gapped bath. 



We have to qualify this statement by taking into account 
the fact that the coupling between electrons and bath al- 
ways also induces an effective interaction between the elec- 
trons. In this way, a given electron becomes coupled indi- 
rectly to the bath of other electrons, such that scattering 
processes will, indeed, lead to a finite decay rate even for 
those low-lying excited levels. This precludes any rigorous 
proof demonstrating the complete absence of decay and 
dephasing for low-lying levels even in such a rather simple 
situation, in spite of the assumption of an excitation gap 
in the bath spectrum. However, the consequences of the 
effective interaction can be distinguished easily from the 
finite dephasing rate that would be predicted by looking at 
Si alone. If the coupling to the bath is of strength g, then 
the latter rate would go as g 2 . In contrast, the effective 
interaction (obtained after integrating out the bath) will 
itself be of strength g 2 , such that the resulting relaxation 
rate (due to coupling of an electron to the density fluctu- 
ations of other electrons) will be of fourth order in g. In 
addition, of course, the rate will depend strongly on the 
distance to the Fermi surface, vanishing when the Fermi 
energy is approached. 
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VIII. CONCLUSIONS 



Appendix A: SOME QUANTITIES FOR THE 
DAMPED OSCILLATOR 



We have tried to demonstrate that the real part Sr of 
the Feynman- Vernon influence action cannot be neglected 
in an analysis of dephasing and decay near the ground state 
of a system, in spite of the fact that, in the simplest sit- 
uations (involving highly excited states and the semiclas- 
sical analysis), it is only the imaginary part Si which en- 
ters the dephasing rate. To this end, we have discussed 
how Sr and Si may cancel each other's effects not only 
in lowest-order perturbation theory but also to all orders, 
by examining exactly solvable linear quantum dissipative 
systems. In general, the cancellation is only found after 
proper integration over the fluctuations away from semi- 
classical trajectories, taking into account the action of the 
unperturbed system and its initial density matrix. Further- 
more, we have pointed out an essential difference between 
the damped oscillator and the free particle with respect to 
this issue. We have argued that the insight obtained in the 
case of the damped oscillator is also applicable to nonlin- 
ear systems and important for discussions of dephasing in 
systems of disordered degenerate fermions. 

In summary, it may be possible to discuss dephasing and 
decay without considering Sr either if the noise is nearly 
classical (external nonequilibrium radiation or bath in the 
high-temperature limit) or if the system itself is in a highly 
excited state. Otherwise, dephasing rates obtained solely 
from Si are bound to come out finite at zero temperature 
in most cases, even when a simpler Golden Rule calcula- 
tion gives vanishing results. Judging from the examples 
discussed above, this is probably not because such a non- 
perturbative procedure goes beyond the Golden Rule ap- 
proximation, but because it neglects some physics already 
contained even within this approximation. 

Still, although the primary message of this paper is that 
special care has to be taken in extracting dephasing rates 
from calculations using the influence functional, we empha- 
size at the same time that there is no general proof showing 
the impossibility of zero-temperature dephasing near the 
ground state. There cannot be such a proof, since there 
is evidently at least one counter-example (the Caldeira- 
Leggett-4H«|dcl of quantum Brownian motion and similar 
modelsEJfeaEa) , and, furthermore, there is no generally ap- 
plicable definition of "dephasing" that is useful under all 
conceivable circumstances. 
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In this appendix, we list, for purposes of reference, the 
quantities relevant to our discussion of the time-evolution 
of the quantum damped harmonic oscillator. These can 
also be found in Ref. [18| (in a slightly different notation). 

The quantity Coo arises in evaluating Si along a pair 
of semiclassical paths, by inserting the solution r(-) of Eq. 
( p2|) for the boundary conditions given by r t , r$ and de- 
termining the coefficient of rg in the result (see Eq. (25)). 
Necessarily, Coo contains only the symmetric part of the 
bath correlator, since Si depends only on that: 



Coo = ~ (sin(wt)) 2 J dti J dt 2 exp + t 2 fj 

x sm(u>{t - ii)) ({ifyi), F(t 2 )}) sin(w(i - t 2 )) . (Al) 



Here Cj = y/uj^ — (t/2) 2 is the renormalized frequency 
of the (underdamped, 7 < 2ui ) oscillator. 

The other quantities needed for calculating (q 2 (t)) arise 
from the evaluation of ReS c i (Eq. (pi|)). In contrast to 
Coo, which has to be evaluated numerically, they can be 
given in closed form: 



771(7/2 + ui cot(u)t)) 



(A2) 

sin(<W) (A3) 

The exponential increase of L m cancels that of Coo when 
calculating the width in Eq. (|29|). However, if the damping 
rate 7 is set to zero, Coo still grows beyond all bounds while 
L t o remains bounded. 



Appendix B: DENSITY MATRIX PROPAGATOR 
FROM THE WIGNER DENSITY EVOLUTION 



Using the classical Langevin equation (|3Cj), the kernel J 
which relates the reduced density matrix of a linear damped 
system at time t to that at time can be found in a way 
which is physically more transparent than the correspond- 
ing derivation using path-integrals (see Eqs. (|]), (||), (|l§|), 
( p4|) and (25)). One first solves for the time-evolution of 
the classical phase space density (i.e. the Wigner density) 
under the action of friction and the Gaussian random force 
F(-). Starting from a 5 peak located in phase space, the 
phase space density evolves into a two-dimensional Gaus- 
sian distribution, whose covariances are related to the cor- 
relator of the force. This gives the propagator J w of the 
Wigner density, which only needs to be Fourier transformed 
with respect to the momenta in order to obtain the den- 
sity matrix propagator J, expressed via center-of-mass and 
difference coordinates R and r: 



/ 



J(Rt,r t \R ,r ;t) = 
dp t dp e^ rt -^J w (R t7 p t \Ro,Po;t). (Bl) 
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Here we show how the propagator J for the density ma- 
trix of a free damped particle subject to the Ohmic bath 
may be obtained in this way. Everything works the same 
for the damped oscillator, only the resulting expressions 
are slightly more lengthy. 

The propagator J w (Rt,pt\Ro,Po',t) of the Wigner den- 
sity is found by solving the classical equations of motion 
for R and p for a given initial condition (Ro,po), taking 
into account friction and the action of the force F: 



^ = - 1P + F(t) 



dp 

m 

dR 
~dt 



P_ 

m 



This yields the solutions 



Pt 
R, 



Poe 7 * + C P 

i?o + ^(l-e-^) 
V 



(B2) 
(B3) 



(B4) 
(B5) 



where £ p and £r are given as integrals over the force F(-): 



& = f dse-^-^Fis) 
Jo 

£r = - f ds(l-e-^- s >)F( S ), 
V Jo 



(B6) 
(B7) 



Since F(-) is a Gaussian random process, £ p and £r are 
Gaussian random variables as well. Therefore, the phase 
space density evolving out of S(R — Ro)S(p — po) is a two- 
dimensional Gaussian distribution in phase space (R,p): 



J w \Rt,Pt\Ro,Po;t) = (6{Rt -Rt- £ R )6(pt -pt- &)} ■ 

(B8) 

The average values p t and R t may be read off from eqs. 
(B4) and (B5). J w has to be Fourier transformed with 



respect to pt and po in order to arrive at the density matrix 
propagator J(i? t , r t \ Ro, ro; t) (see Eq. (Bl)). In order to 
do this, we express as a Gaussian density in terms of 
Pt, Po, for fixed R t ,Ro: 



Making use of (B4) and (B5), the average values and the 
deviations of pt and pq are found to be given by: 



Pt 
Po 
Spt 
Spo 
A 



X(Rt — Ro) 
e*\(Rt - Ro) 

V 



e7 t 



1 



(B12) 
(B13) 
(B14) 
(B15) 

(B16) 



po and pt are the momenta at time t and which the 
particle must have if it is to go from R to R t in time t, 
provided no fluctuating force acts. Spt and Spo are the de- 
viations from these values needed to compensate the effects 
ofF(-). 

Quadratic completion in the exponent immediately 
yields the result of the Fourier integration over pt and po, 
which is the desired density matrix propagator J: 



J(Rt,r t \Ro,r ;t) oc 



exp 



i{R t - Ro)X(r t - e^r ) - ~ {ir t 8p t - r 5p f) 



B17) 



This reproduces the result given in Ref. [is] (or Ref. 0, 
App. E). The prefactor can be determined from the nor- 
malization condition, Eq. (20), and only depends on the 
time t. The term in angular brackets still has to be aver- 
aged over the force F(-). The resulting real part of the ex- 
ponent equals — Si evaluated along the semiclassical paths 
(compare the general structure given in Eq. (|25[)). In terms 
of F(-), Spt and Spo read explicitly: 



S P t = 
Spo = 



ds ■ 



ds ■ 



-F(s) 



aft 



err* 



-F(s) 



(B18) 
(B19) 



J oc exp 



with 



—5P t K- 1 6P 
2 



SP = 



Pi 
Po 



and the covariance matrix 



K = 



{5pi) 
{SptSpo) 



Pt 

- Po 



(SptSp ) 
(5p 2 o) 



(B9) 



(BIO) 



(Bll) 



Note that 



Spt ~8po = / dsF(s) 



f 

Jo 



(B20) 



The avera ges (6p$), (Spo) and (Sp t Sp ) to be evaluated 
in Eq. (B17) contain time-integrals over the force correla- 
tor (F(si)F(s2)) . At zero temperature, these integrals lead 
to terms growing logarithmically in time, which are char- 
acteristic for the free particle coupled to an Ohmic bath 
(compare the discussion in the main text, sec. VI). 
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